Analysis on the effect of technical fluctuations on laser lineshape 
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We analyze theoretically the effect of technical fluctuations on laser lineshape in terms of statistics 
of amplitude and phase noise and their respective bandwidths. While the phase noise tends to 
broaden the linewidth as the magnitude of the noise increases, the amplitude noise brings out an 
additional structure with its spectral density reflecting the magnitude of the noise. The effect of 
possible coupling between those two noises is also discussed. 

PACS numbers: 



I. INTRODUCTION 

Lineshape is one of the fundamental properties which 
tells about how monochromatic a laser is. It also provides 
a way to investigate the light emission mechanisms for 
various types of lasers [lj . However the technical fluctua- 
tions inherent in a laser and its environment prevent one 
from observing spectral structure determined by intrinsic 
dynamics only. In fact, the intrinsic limit of linewidth set 
by quantum noise is much smaller than we usually ob- 
serve in laboratory. The sources of such noise is usually 
thermal or mechanical and they appear as, for example, 
cavity drift, fluctuation of population inversion and insta- 
bility pump field. We may classify the aspects of those 
fluctuations into phase and amplitude noise with vari- 
ous bandwidth. We analyze the effect of technical fluc- 
tuations on lineshape in terms of amplitude and phase 
in case of sinusoidal modulation and Gaussian distribu- 
tion. Their respective bandwidth and possible coupling 
between those two Q are also included. 



II. PHASE NOISE 

Let us write the time variation of the electric field, the 
spectrum of which we want to measure, as 



E(t) = E {t)e i ^ ot+ ^ 



(1) 



where Eo(t) is a slowly varying envelope and <p(t) is a 
randomly fluctuating phase. 

Firstly consider the case in which the amplitude E$ (t) 
is constant Eq and 4*{t) undergoes a random walk process 
which imposes Gaussian statistics on cj)(t). The effective 
frequency is 



(2) 



which represents that the time derivative of (f>(t) makes 
jitter around the carrier frequency. According to Wiener- 
Khinchin theorem, the spectral lineshape is given by 



Fourier transform of the first-order correlation function 



g a ){T) _ (E*(t)E(t + T)) 



(3) 



where ( ) denotes a time average. Using Eq. (1), 

(E*{t)E(t + r)) = |Eb| a e toor (e i tf +T ^*'> ,ft ') . (4) 

For a normally distributed random variable x with its 
probability density function (PDF) 



P[x] = 



^_ e -(x-x) 2 /2<r 2 



2n<r 



(5) 



where x is the mean of x and a 1 is the variance, 
the characteristic function (e lux ) is calculated to be 
e ixu e -\o*^ Taking x = J* +T <j>(t')dt' and lu = 1, 



{E*{t)E{t + T)) = \E Q \ z e^> T e 
where x — and a^. is given by 



t + T rt + T 

dt' / dt"^(t')4>(t")) 
t Jt 



(6) 



(7) 



The calculation of cr^ depends on the bandwidth of the 
frequency jittering. 

It is readily calculated if the spectrum of <j>(t) is flat, 
that is, (<j)(t)<j)(t + t)) is proportional to <5(r). The con- 
stant of proportionality is found from 



4> 2 (t)p[<p(t)}d<p(t) 

1 



V2^AR 2 
(AR) 2 , 



■(2Ai? )' 



(8) 



where AR denotes the magnitude of jittering and the 
PDF of 0(f) 
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FIG. 1: (a) The correlation (<p(t)(j>(t + r)) as a function of 
r as a model of finite bandwidth Gaussian noise, (b) The 
value of integration of Eq. (10) is the area surrounded by 
t" =t'± V2,t' = r and t" = r 



was used. Then 

2 



dt' / dt" {4>{t')<j)(t")) 
o Jo 

= (AR) 2 f dt' f dt"5(t' - 1") 
Jo Jo 

= V2(AR) 2 \r\ , 



(10) 



where, in the first line, t is replaced with zero regarding 
the noise being stationary. Finally we came to obtain 



■73( A «) 2 I-I 



(11) 



and its Fourier transform gives a Lorentzian line shape 
with FWHM of (AR) 2 . 

If the spectrum of <j)(t) has a finite bandwidth, a 2 de- 
pends on the correlation time t c of (f>(t) . We can model it 
by a square-shaped temporal correlation as in Fig. 1(a). 
In terms of frequency this noise is the white noise filtered 
by a sine- function low-pass filter of widths t~ . The in- 
tegration in Eq. (10) then takes the value depending on 
r(Figl.(b)) 



(ARY 



T 



2t 2 



if t < V2t c 
if t > y/2t c 



(12) 



Thus g^(r) is Gaussian up to r = \/2t c and thereafter 
exponentially decaying function. Of course it is continu- 
ous at t = \[2t c having the common value e~( 2Ai? ) '<=. If 
t c is longer than (AR) -1 , g^\r) is practically Gaussian. 
We have Gaussian lineshape with its linewidth ~ AR in 
that case. In the opposite limit where t c goes to zero, the 
Lorentzian lineshape due to the white noise is recovered. 

Another possible case of frequency jittering is slow 
modulation of carrier frequency as it usually happens in a 
laser with its cavity slowly drifting around the resonance 
frequency. PDF of Eq. (9) is no more valid in such cases. 
Rather we start from 



<j)(t) = AFcos(flt) 



and accordingly 



AF 

IT 



• sin fit = 4>q sin fit , 



(13) 



(14) 



where AF is the amplitude of modulation and fl is the 
slow frequency. For this simple harmonic oscillation, 
PDF of <j)(t) is given by 



PW)] 



1 



1 



(15) 



The expectation value of e zx will be calculated using this 
PDF where 



t + T 



4>(t')dt' 



t + T 



AFcoaflt'dt' 



6o[sin Cl(t + t) — sin fit] 
bo (sin fit cos f2r + cos f2i sin f2r — sin fit) 
f2r s 



4> x (—2 cos — ) ± sin fir 



(16) 



By inspecting the PDF in Eq. (15), we can recognize 
that 4>(t) spends most of its time near <fi(t) ~ (j)Q thereby 
neglect the second term in Eq. (16). Then 



00 



-2id> cos 



' c le 



-2id) cos 



00 T - (j) 2 



Jo 20o cos 



(17) 



where the integration in the last line involving Jo , Bessel 
function of the first kind is performed in Ref. . There- 
fore 



. . 9 fir 
J [ 20o cos z — 



(18) 



The graph of Jo(20o cos 2 rir/2) is given in Fig. 2(a) 
and (c) for different value of magnitude of modulation 
AF. It is infinite pulse train with the repetition rate 
(f2/27r) _1 . Its Fourier transform constitutes frequency 
comb whose width of envelope is determined by the in- 
verse time duration of the pulse, T. (Fig.l (b),(d)) T 
gets shorter as we increase AF because the argument of 
Bessel function changes by larger amount for the same 
change of r. This is reasonable result since the more 
harsh we swing the carrier frequency, the wider the spec- 
trum should be. f2 determines the repetition rate, i.e. 
density (degree of fine-tooth) of the comb. 

If another frequency components other than f2 is added 
in Eq. (12), the fine spectral structure is easily destroyed 
by the complexity of argument in Bessel function in Eq. 
(17). The additional modulation is linearly added to 
x and finally to the argument of the Bessel function. 
This complicated argument brings about reduction of the 
pulse height of Jo by only a small number of such super- 
position. Thus, in the spectrum, the spectral density of 
each tooth in the comb gets smaller and interval between 
adjacent teeth gets narrower. 
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FIG. 2: (a) The factor J o (20 o cos 2 ^) appearing in Eq. (18) 
as a function of r. (b) Fourier transform of (a), (c) The factor 
Jo (2<^)o cos 2 Qf) for AF 10 times larger than that in (a), (d) 
Fourier transform of (c). 



III. INTENSITY NOISE 

Next let us include the effect of the amplitude fluctua- 
tion. The criterion between amplitude and phase, in this 
analysis, rests with their direct appearance in intensity. 



(a) 




FIG. 3: (a) Sidebands due to the amplitude modulation, (b) 
If the modulation has a finite band, the resultant spectrum 
has low-lying wings corresponding to the noise band. In this 
figure, Gaussian noise is presented. 



Consider the intensity modulated like 

\E Q (t)\ 2 = I (l + M cos fit) 



(19) 



where the modulation depth M is usually much less than 
1. The corresponding field amplitude can be written as 

E (t) = v^o(l + aie mt + a-ie"** + a 2 e l2nt ■■■), (20) 

where 1 » a\, <z_i and 01, a_i » 02 ■ ■ ■ for M <C 1. This 
implies that the spectrum is given with sidebands, which 
are symmetrically apart from the carrier frequency by 
integer multiples of Q. (Fig. 3(a)) Extending this idea to 
the noise with finite band we can imagine superposition 
of many sidebands comprised of all frequency component 
within the band. The resultant spectrum contains the 
low-lying wing structure near the carrier frequency. One 
example where the Gaussian noise is applied is depicted 
in Fig. 3 (b). Note that the Lorcntzian is smoothed out 
by Gaussian profile which brings about the considerable 
deviation from a Lorentzian lineshape. 
The explicit calculation starts with 

(E*{t)E{t + r)) =e WoT (E*(t)E (t + T)e l It t+r ^ t ' )dt ') . 

(21) 

Since 



E*(t)E (t + r) 



h(l + \a 1+ \ 2 e+^ 



each sideband make x to shift by 0, O, 2f2 ■ 
while ol remains the same. Thus 



\a\- 
■), (22) 
• respectively 



9 {1) (r) = 
so that 



e iW()T e -^(A-R) T 

x (l + \a 1+ \ 2 e +inT 



ai_ 



+ 



(23) 



(uj - io r + (ARyVz) 2 



+ 



\ a l+\ 



(a)-^o-«) 2 + (Ai? 2 /V2) 2 



(24) 
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The same result can be understood in a different way. 
The spectral amplitude £(u>) is the Fourier transform of 
the product of two functions: E (t) and e *M+<K*)]. Thus 
the overall spectrum is given by convolution of each spec- 
trum of the functions. If the bandwidth of the amplitude 
noise is Gaussian, 



5(u' - ujo) + ace 



(oj'-oj ) 2 /A g 2 



(lj'-lj- tu y + A| 



■ duj' , (25) 



where A t, and Aq are spectral widths of the Lorentzian 
and Gaussian respectively and aa is relative spectral am- 
plitude of Gaussian noise. 

If the bandwidth of amplitude fluctuation is too broad, 
the low-lying structure may not be easily recognized. The 
manifestation of amplitude noise can then be confirmed 
in the experiment from the intensity correlation function 
g( 2 \r) since 

.9 (2) (0) 

= ((1 + J m{fl) cos fit dfl)(l + J m{fl') cos fl'tdfl')) 
= 1 + ^Jdflj dflm{fl)m(fl')S(fl - fl') 

= 1 + \J dnm2 W ' ( 26 ) 
where m(u>) is modulation density. The peak near r = 



definitely reveal the effect of amplitude fluctuation. 

There might be possible coupling between amplitude 
and phase as is well known in semiconductor lasers. If 
4>{t) is affected by the intensity modulation as a quadra- 
ture 4>{t) = 28 sin fit where 8 designates the magnitude 
of coupling, 

E(t) ~ (1 + 2cn cos fit) cos (ut + 25 sin fit) 

= ^e wt + (ai + <5)e J ^+°) + (ai - ^e"*^-^ 

(27) 



+c.c. 



Hence this type of coupling brings about the asymmetry 
in the sidebands. 



IV. CONCLUSION 

We investigated the effect of technical fluctuations on 
laser lineshape. The Gaussian noise in frequency make 
the lineshape a Lorentzian or a Gaussian depending on 
the correlation time of the noise. Slow swing of frequency 
results in a lineshape with many sidebands by the period- 
icity involved. The amplitude noise imposes the low- lying 
wing structure in the spectrum. The bandwidth of the 
noise only determines the width of the additional struc- 
ture while the coupling between the amplitude and phase 
might lead to an asymmetry in the structure. 
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